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Recently it has been shown that multicomponent spin-orbit-coupled fermions in one-dimensional 
optical lattices can be viewed as spinless fermions moving in two-dimensional synthetic lattices with 
synthetic magnetic flux. The quantum Hall edge states in these systems have been observed in 
recent experiments. In this paper we study the effect of an attractive Hubbard interaction. Since 
the Hubbard interaction is long-range in the synthetic dimension, it is able to efficiently induce 
Cooper pairing between the counterpropagating chiral edge states. The topological class of the 
resultant one-dimensional superfluid is determined by the parity (even/odd) of the Chern number 
in the two-dimensional synthetic lattice. We also show the presence of a chiral symmetry in our 
model, which implies Z classification and the robustness of multiple zero modes when this symmetry 
is unbroken. 

PACS numbers: 73.43.-f, 71.10.Hf, 74.20.Mn, 71.45.-d, 71.70.-d 


Topological superconductors and topological superfluids hosting Majorana zero modes [1-7] have been among the 
central themes of both condensed matter and cold atom physics recently. (“Topological superconductor” refers to 
charged particles, while “topological superfluid” refers to neutral particles, otherwise their physics is essentially the 
same. The result of our paper is equally applicable to topological superconductors and topological superfluids.) Apart 
from being novel phases of matter, they have potential applications in quantum computation [8, 9]. It is therefore 
highly desirable to search for various routes towards realization of topological superconductivity/superfluidity and 
Majorana zero modes. There have been several proposals to realize them in either condensed matter [10-14, 14-19] 
or cold atom systems [20-23]. The latter systems have the advantage of high controllability. Experimental study on 
topological superconductors and Majorana zero modes is also extremely active [24-28]. 

In this paper we study the Cooper pairing between chiral edge modes of quantum Hall strips as a possible route 
toward one-dimensional (ID) topological superconductors and topological superfluids. This is stimulated by the idea 
of “synthetic dimension” [29, 30] emerging from optical lattice with atoms with large spin[31-42]. In this visualization, 
internal degrees of freedom (“spin”) form an additional spatial dimension. This picture is especially convenient when 
the internal states are coupled sequentially, as can be readily done by two Raman beams [29]. Synthetic magnetic 
flux naturally exists inside the two-dimensional (one physical dimension plus one synthetic dimension) lattice, and 
the quantum Hall states (Chern insulator) can be simulated. The chiral edge modes have been observed in recent 
experiments [43, 44]. The two counterpropagating chiral modes are separated in the synthetic dimension, therefore 
they are immune to backscattering if the scatters are short-range. 

The interaction effect in this system is an important issue[29, 45-47], both experimentally and theoretically. With 
the motivation of realizing topological superfluidity, in this paper we study the effect of an attractive Hubbard 
interaction, which is long-range in the synthetic dimension, though short-range in the physical dimension. As a 
consequence, the two counterpropagating chiral edge modes can be Cooper-paired efficiently. We find that topological 
superfluidity naturally arises. We also study the existence of topological Majorana zero modes. Interestingly, multiple 
zero modes are stable if a chiral symmetry of our model is unbroken. 


Results 

The model and main picture. The system in consideration is illustrated in Fig.l. The one-dimensional optical 
lattice extends in the x direction, and the Raman-induced hopping couples M = 2F + 1 spin states at each site x in 
a sequential manner (F = 5/2 is shown in Fig.l). Therefore, the system acquires a “synthetic dimension” [29] . The 
Hamiltonian is = Hq -j- Hi, in which the free part is 


Ho 
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where Pirn = CtgF,m = Ct^/F{F + 1) — m{m + 1)[29], with Cl depending on the strength of Raman transitions. Gener¬ 
ally, the value of Utrn can be controlled by tuning Raman beams. Our results will be insensitive to details of ^2^. The 
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FIG. 1: Sketch of the system. The spin states are coupled sequentially by Raman-induced hopping, which generates a synthetic 
dimension, in addition to the physical dimension x. 


above Hamiltonian is readily realizable in experiment [43, 44]. The presence of indicates that there is a flux 7 

in each plaquette, which is responsible for the emergence of the chiral edge states in this model(see Fig.(l)). After a 
gauge transformation Cn,m Hq becomes 

^0 = + H.C.) 


^n,m^n,m5 (2) 

n,m 

in which the hopping gains a spin-dependent phase factor. 

In this paper we take the simple yet realistic SU(M)-invariant Hubbard interaction: 

= (3) 

n 

where Nn = Xlm m^ri,m and U < 0. This interaction is apparently long-range in the synthetic dimension, thus it is 
quite capable of pairing the counterpropagating modes at the opposite edges (near m = F and m = —F, respectively) 
in the synthetic dimension. 

Before proceeding to a quantitative study of iJ/, we would like to discuss the physical picture of possible topological 
superfluidity in this model. Since the hopping lacks translational symmetry along the synthetic dimension, the bulk 
Chern number as an integral of Berry curvature in the two-dimensional Brillouin zone cannot be defined, however, its 
manifestation as the number of chiral edge modes is well-defined. Suppose that the “bulk Chern number” (7 = 1, i.e. 
there is a single pair of chiral edge modes in the bulk gap [see Fig. 2, in which the purple dotted lines intersect with 
the chiral modes]. If there is a small Cooper pairing between these two edge modes, the system is a one-dimensional 
topological superfluid. This can be inferred using Kitaev’s Z 2 topological invariant [2] , which essentially counts the 
parity (even/odd) of the number of Fermi points within [0,7r/a] in the absence of pairing. In the case (7 = 2, i.e. 
there are two pairs of chiral edge modes, as illustrated by the curves intersecting with the blue dotted line in Fig. (2b), 
the superfluid resulting from pairing the edge states is topologically trivial. Provided that the pairing is small, the 
superfluid (or superconductor) is nontrivial (trivial) when the bulk Chern number is odd(even), namely, 

(-1)*^ = (-1)^, (4) 


where z/ = 0,1 (mod 2) is the Z2 topological number of ID superconductor or superfluid [2]. In the rest part of this 
paper we shall present a quantitative study of the picture outlined above. 


Cooper pairing in self-consistent mean-field. At the mean-field level the Hubbard interaction can be decomposed 
as 


Hi = 
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where we have defined An,mm' = U < Cnm'Cnm >, which satisfies An,mm' = —An,m'm because of Fermi statistics. 
The mean field BdG Hamiltonian becomes 


Hmf = Ho + 1 
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FIG. 2: The spectrum of free Hamiltonian Hq for (a) F — 5/2, t — 1,0. — 0.1, 7 = tt/S. The dotted line is located at /i = —1.6 
; and (b) F — 7/2, t = 1, O — 0.1, 7 = 7 r/ 3 . The purple dotted line is located at /xi = —1.8, and the blue dotted line is located 
at [12 = —0.95. 




FIG. 3: The self-consistent solutions of Gooper pairing Amm' as functions of U for (a) = —1.8; (b) (12 = —0.95. The 

parameters are t = 1, Q = 0.1, and 7 = 7 r/ 3 , which are the same as used in Fig. (2b). 


Since the basic physics is the pairing of chiral edge modes with opposite momenta, it is natural to consider Cooper 
pairing with zero total momentum, namely that An,mm' = ^mm' is independent of n. We determine these Amm' 
self-consistently. Whenever there are several sets of self-consistent solutions of {Amm'}^ we compare their mean-field 
energies and pick up the ground state. This calculation can be carried out for all possible values of F. Hereafter we 
take F = 7/2 as an example. 

The pairings as functions of Hubbard U are shown in Fig. (3) for two values of chemical potential fi. In Fig. (3a) 
we take /ii = —1.8, which corresponds to bulk Chern number C = 1 [see the purple dotted line in Fig. ( 2 b)]. For 
small \U\, as \U\ increases, Amm' grows exponentially. For U not too large, A^ 12,-112 dominates other pairings, 
which is consistent with the picture of pairing between chiral edge modes. At a critical \U\ slightly below 2.5, the 
system undergoes a first-order transition to a phase in which Ai/ 2 ,- 1/27 which should be regarded as a “bulk pairing”, 
becomes comparable to the edge pairing. Thus the edge-pairing picture becomes inaccurate at large \U\. In Fig. (3b) 
a different chemical potential 112 = —0.95 is taken, and the edge-state pairings again dominate at small \U\. 

To illustrate the robustness of the present picture, we study a set of quite different parameters, shown in Fig. (4). 
The behavior of pairings is qualitatively the same as that for the previous parameters. 

A few remarks are in order. First, although we have taken the strength of hopping along the synthetic dimension 
as OgF,m, we have also checked that the result is qualitatively the same when it is m-independent. Second, when the 
2D bulk is metallic, topological superfluidity can still emerge, though there is no clear criteria using Chern number. 
We shall not focus on details about this. 

Majorana zero modes. The hallmark of ID topological superconductor or superfluid is the emergence of topological 
zero modes localized near the two ends of an open chain. 

We have solved the BdG mean-field Hamiltonian Eq.( 6 ) for the wavefunction ip{x,m) = [u{x,m),v{x,m)Y, in 
which u{x,m) and v{x,m) denotes the particle and hole component respectively. Below we present our solutions in 
a chain with sharp boundary, for the parameters t = 1, O = 0.1, 7 = 7 r/ 3 , at both /ii = —1.8 and /i 2 = —0.95. The 
Cooper pairings are taken to be the mean-field values at U = —1.3, which we obtained in the previous section. The 
case 111 = —1.8 is shown in Fig. 5. There is one zero mode localized at each end of the open chain, and a tiny finite-size 
coupling mixes them slightly, though the energy splitting due to finite-size effect is too small to be discernable. The 
existence of a single Majorana zero modes at each end of an open chain is consistent with Eq.(4), the bulk Chern 
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FIG. 4: (a) The spectrum without interaction for t = 1 , Q = 0.025, and 7 = tt/G. The dashed line marks fi = —1.94. (b) The 

Cooper pairings as functions of U. 



FIG. 5: (a) The wavefunction of a zero mode in an open chain with length N = 150. The parameters are t = l,Q = 0 . 1,7 = 7 r /3 
and 11 = —1.8. The Cooper pairings are taken as Ayi 2 ,- 7/2 — 0.09, A 7 / 2 ,- 5/2 = —0.019, A 5 / 2 ,- 5/2 = 0.03, A 5 / 2 ,- 3/2 = 
—0.007, A3/2,-3/2 = 0.01, A3/2,-1/2 = —0.002, Ai/ 2 , -1/2 = 0 . 002 , which are the mean-field pairing obtained at \U\ = 1.3 
( see also Fig. (3a) ). Other A^^/’s are much smaller and thus neglected. The two zero modes have the same profile of 
\'ip(x)\^ = |'0(x,m)|^, thus only one is shown here. The inset shows a few energies near E = 0. (b) The zero mode solution 

with m resolution. 


number being (7=1 (odd number) at jxi = —1.8. The superfiuid is topologically nontrivial in this case. 

As a comparison, we also present the zero mode solutions at /i 2 = —0.95, for which the free Hamiltonian have two 
pairs of chiral edge modes ((7 = 2), indicating that the superfluid phase at small \U\ should be Z 2 topologically trivial 
(see Eq.(4)). In the numerical calculation with open boundary condition, we find two Majorana zero modes at each 
end (see Fig. 6), which means that the superfluid is Z 2 trivial. Therefore we see again that (—1)^ determines the Z 2 
topological classification of the ID superfluid in synthetic dimension. 

One may wonder why there is no hybridization between the two zero modes, which may open a gap for them. We 



FIG. 6 : The wavefunction of zero modes in an open chain with length N = 200. The parameters are t = l,f7 = 0.1, 7 = 7 r /3 
and = —0.95, which is the same as Fig. 5 except for the chemical potential. The pairings are Aj/ 2,-712 — 0.067, A 7 / 2 ,- 5/2 = 
—0.005, A 5 / 2 ,_ 5/2 = 0.052, A 3 / 2 ,_ 3/2 = 0.028, Ai/ 2 ,_i /2 = 0.009 (the mean-field pairing at (7 = —1.3), and other pairing terms 
are much smaller and thus neglected. The spatial profile of |'0(x)|^ of two of the four zero modes are shown in (a) and (b), 
while the profile with m resolution is shown in (c) and (d). The other two zero modes with the same profiles are not shown 
repeatedly. The inset of (a) shows several energies close to E = 0. 
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FIG. 7: Majorana zero modes in the harmonic trap. The parameters are t = 1, Q = 0.1, and 7 = vr/S, which are the same 
as used in Fig. (2b). The harmonic trap is V{x) = 0.00025(x/a)^, where x = 0 is the center of the chain with size N = 150. 
The two zero modes have the same profile of |'0(x, ?n)|^, thus we only show one. The inset of (a) shows several energies close 
to = 0 . 


shall explain the reason as follows. In fact, the BdG Hamiltonian has a time-reversal symmetry and a particle-hole 
symmetry, which can be combined into a chiral symmetry [49]. If the Cooper pairing are real, then we can check 

that the Hamiltonian satisfies 


in which 


CHuF{k)C-^ = -HuF{k) 

( 

—i 

C = 

i 

V* / 


( 7 ) 


(8) 


This matrix is written in the BdG basis of {ck^F^ * * * ^Ck,-F^ F’ ''' ’ ^-/c -f)* these symmetries, the BdG 

Hamiltonian can be classified as BDI[48, 49], whose classification in ID is Z. The Z 2 topologically trivial phase is 
nontrivial according to the Z classification of BDI class, which is the reason why zero modes appear at the edge of 
Z 2 trivial states. We have checked that, if we break the symmetries, e.g. by giving a phase factor to A 7 / 2 ,- 7/2 (with 
other Amm'S unchanged), then these zero modes will be shifted to nonzero energies. 

According to Fidkowski and Kitaev’s work[52], in the presence of interaction, the classification of BDI-class topo¬ 
logical superconductors in ID is Zg instead of Z. Because of the flexible tunability to topological superconductors 
with large topological number (say 8) in our system, the Zg classification can be tested experimentally. If we tune 
the bulk Chern number of our system to 8, the eight nominally-zero-modes will be shifted away from zero energy due 
to the (beyond-mean-field) interaction effects of the Hubbard term. If observed, this will be an experimental test of 
Fidkowski and Kitaev’s Zg classification. 

To make a closer connection to experiment, we also study the existence of Majorana zero modes in a system with soft 
boundary created by a harmonic trap V(x). In the presence of V(x), the chemical potential becomes ju(x) = jXo — V{x). 
We take jno = —1.7 and V{x) = 0.00025(x/a)^, such that the center of the system is topologically nontrivial. Since 
fi{x) is not constant, Amm' should also be x-dependent. To incorporate this effect, we numerically calculate the 
functions A^^/(/i) at = —1.3, which is then used to produce the mean-field BdG Hamiltonian in harmonic trap. In 
the solution to this BdG Hamiltonian, the zero modes can be clearly seen, as shown in Fig. (7), though the quantitative 
details are different from the case of hard boundary. 

Conclusions and Discussions. 

We have studied the pairing between counterpropagating chiral edge modes in the quantum Hall strip in synthetic 
dimension. This picture has several merits. Creation of magnetic flux in the synthetic dimension by Raman beans is 
easier than in physical dimensions. The spatial separation of left and right moving chiral edge states in the synthetic 
dimension effectively prevents the backscattering between them, which implies their robustness. Meanwhile, the 
Hubbard interaction is not suppressed by this spatial separation: It is infinite-range in the synthetic dimension, 
therefore, it can pair the two edge modes quite efficiently. As we have shown, the resultant states are topological 
superfluid carrying Majorana zero modes. If the chiral symmetry of our model is unbroken, the classification is Z and 
multiple zero modes are stable; on the other hand, if this symmetry is broken, the classification is Z 2 . 
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Finally, we remark that quantum fluctuations of the phase factor of pairing in ID is generally strong. One can put 
the ID system in proximity to a 3D supefluid to suppress these fluctuations [20]. Moreover, it has been shown[50, 51] 
that long-range superconducting order is not a necessary condition for the existence of Majorana zero modes. The 
zero modes persist even when the long-range superconducting order is replaced by algebraic order (i.e. the correlations 
of pairings decay by power-law). In our system this conclusion applies. 

Methods 

Mean-field calculations. The mean-field calculation is carried out by the standard procedure of decomposing the 
Hubbard interaction as fermion bilinear terms, leading to Eq.(6). The Cooper pairing is calculated from Eq.(6) in a 
self-consistent manner. All self-consistent solutions for the Cooper pairing are obtained. In the case that there are 
more than one self-consistent solutions, the one with lowest mean-field energy is selected. 
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